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VACUUM SOLUTIONS WITH NON-ALIGNED 

POLARIZATION 



The purpose of this chapter is to describe the known exact solutions for 
colliding gravitational waves in which the polarization of the approaching 
waves is not aligned. The first solution of this type was obtained by Nutku 
and Halil (1977). A further generalization of this solution was attempted 
by Halil (1979), but this has subsequently been found to be incorrect and 
so will not be considered here. 



13.1 The Nutku-Halil solution 

The Khan-Penrose solution discussed in Chapter 3 describes the collision 
of impulsive waves whose polarization vectors are aligned. Nutku and 
Halil (1977) have generalized this solution to give one which describes the 
collision of impulsive gravitational waves with non-colinear polarization. 
This is the most simple solution of this type. 

In terms of the Szekeres line element (6.20), the metric functions of 
this solution can be written in the form 

e- u = l-u 2 -v 2 



e -v ^ /(I -£)(!-£) 



sinh W = 



e -M = 



(1 + E)(1 + E) 
i{E-E!) ( 13J ) 



1 — EE 

1 — EE 



\J\ — u 2 \/l — v 2 Vl — u 2 — V 2 



where 



E = e lc Wl - v 2 + e ip v^l - u 2 (13.2) 



and a and (3 are constants such that (a — (3) is the angle between the 
polarization vectors of the approaching waves. This solution can be seen 
to reduce to the Khan-Penrose solution (3.9) when a = (3 = 0. It may 
also be observed that the functions / and g are given by / = \ — u 2 , and 



g=\-v 2 . 
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This solution has been analysed in great detail by Chandrasekhar 
and Ferrari (1984). Their approach commences by rewriting the field 
equations in the form of Ernst's equation as described in Chapter 11. It 
can then be seen from (11.20) that E is the associated Ernst function. 
Chandrasekhar and Ferrari were then able to show that the Nutku-Halil 
solution is obtained by the simple choice of Ernst function given by 

E=pt + iqz (13.3) 

where p and q are real constants which satisfy the condition p 2 + q 2 = 1. 
It is always possible to use a rotation of the type (12.9) to put (3 = —a, 
and in this case p = cos a and q = sin a. 

It can easily be shown that the Ernst function Z associated with this 
solution is given by 

1 + t cos a + iz sin a 

A = 

1 — t cos a — iz sin a 

It may be observed that (13.3) is in fact the Ernst potential which 
leads to the Kerr solution for stationary axisymmetric space-times. This 
reduces to the potential for the Schwarzschild solution when q = or 
a = 0. Thus the difference in the polarization of the approaching waves 
in this situation can be seen to be equivalent to the rotational parameter 
in the Kerr solution. 

Taking the Ernst functions as (13.3) or (13.4), it can be shown that 
the boundary conditions described in Section 7.2 can only be satisfied if 
f = 2 ~ ( c i u ) 2 + • • •) an d g = \ — (c2v) 2 + — It is then convenient to 
transform the null coordinates such that / = ^ — -u 2 , and g = \ — v 2 . 
Thus the form (13.2) is obtained uniquely from the Ernst function (13.3). 
It also follows from this that the approaching waves necessarily contain 
impulsive components. 

Expressions for the components of the Weyl tensor ^q, ^2 and ^4 
describing the wave and interaction components have been given by Chan- 
drasekhar and Ferrari (1984), though in a notation which differs slightly 
from that adopted here. However, it can still be clearly seen that this 
solution is a generalization of the Khan-Penrose solution in which the 
approaching impulsive waves have different polarization. Apart from the 
initial impulsive waves on the boundaries u = and v = 0, the interiors 
of regions II and III are flat. 

It can also be seen that the singularity structure of this solution is 
the same as that of the Khan-Penrose solution. There is a curvature 
singularity on the space-like surface u 2 + v 2 = 1 in region IV, and there 



(13.4) 
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are fold singularities in regions II and III which are identical to those 
described in Section 8.3. 

It may also be pointed out that Ernst (1986) has shown that that 
the Nutku-Halil solution can also be obtained using a double-Harrison 
(Backlund) transformation, with the seed metric being the isotropic Kas- 
ner metric. 



13.2 The Panov solution 



Panov (19796) has applied the generation method of Geroch (1971), de- 
scribed in Section 12.3, to obtain a non-colinear solution from a colinear 



one. 



In terms of the metric functions in (6.20), Panov started with a co- 
linear solution U , V , and M Q of (6.22) with W = 0. He then worked 
with the Killing vector £ M = 6% , which is (12.25) with a = 1 and 6 = 0. 
This enabled him to put P2 = 0, and hence u = 0, and also from (12.30), 
Q 2 = F 2 - 1 and Q 3 = 0, where F 2 = e v °~ u °. 

With these expressions, the metric components may be derived from 
(12.23), giving 



U = U 

e -2V = e -2V a ^ CQS 2 a + gin 2 a e 2(V -U o) y + ^2 2a p 2 

sinh W = sin2aP 3 e y ° (cos 2 a + sin 2 a e 2{y °~ Uo) j 
^(cos^ + sin^e 2 ^-^)) 



(13.5) 



e -M = e -M a 



where, from (12.28), P3 must satisfy 

Ps,u = e- u °(V Q , u - U , u ) , P s , v = -e- u °(V Q , v - U , v ) (13.6) 

or, more conveniently, 

Psj = (f + g)V oJ + 1, Ps, g = -(/ + 9)V , g - 1. (13.7) 

It is worth pointing out that, if the initial Killing vector had been 
taken to be £ M = 5%, which is (12.25) with a = and b = 1, an equivalent 
transformation would have been obtained with P 3 = and P2 7^ 0. This 
would have given an identical set of equations, but with V being replaced 
by -V, and P 3 by P 2 . 
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In terms of the Ernst formulation of the problem described in Chap- 
ter 11, the above transformation implies that, if Z a is a real solution of 
(11.18), then a new complex solution is given by 

Z = cos 2 aZ + (1 -t 2 )(l - z 2 )sin 2 aZ' 1 + isin2aP 3 (13.8) 
where Ps(t,z) must satisfy the equations 



Panov then took the initial solution to be the totally general solu- 
tion of Szekeres (1972), which can be considered to be equivalent to the 
general solution given here by (10.16). He has thus determined a gen- 
eral class of non-colinear solutions, and he was able to show that, when 
the appropriate boundary conditions have been applied, the approaching 
waves necessarily have variable polarization. 

It has already been shown that the new solution obtained using a 
Geroch transformation will automatically satisfy the required boundary 
conditions if these are satisfied by the initial solution. In this case, it can 
also be seen that the new and initial solutions have the same singularity 
structure. 

13.3 The Chandrasekhar-Xanthopoulos solution 

Chandrasekhar and Ferrari (1984) have shown that the associated Ernst 
function which is contained in the Nutku-Halil solution can be written 
as E = pt + iqz (13.3), where p 2 + q 2 = 1. Now, it is well known 
(see Chandrasekhar, 1983) that this is the Ernst potential that leads to 
the Kerr solution for stationary axisymmetric space-times. However, as 
described in Section 12.4, Chandrasekhar and Ferrari have also shown 
that, for the colliding plane wave situation, the Ernst function may be 
regarded either (a) as containing the metric functions, or (b) as a potential 
for those functions. 

When the Ernst function (13.3) is considered as containing the metric 
functions, it leads to the Nutku-Halil solution described above. However, 
when the same function is considered as a potential, it then leads to 
the Kerr solution. Thus part of the the Kerr space-time must also be 
considered as a solution of the colliding plane wave equations, though in 
this situation the coordinates have a different interpretation. 

This alternative interpretation of part of the Kerr space-time has 
been described in detail by Chandrasekhar and Xanthopoulos (19866). 



P 3 , t = -(l-z 2 )(\ogZ ), z -z 
P 3 , z = -(l-t 2 )(logZ ), t -t. 



(13.9) 



122 Vacuum solutions with non-aligned polarization 

Essentially it may be considered to have been derived from the Nutku- 
Halil solution using the Neugebauer-Kramer involution described in Sec- 
tion 12.4. 

We start here with the Ernst potential given by E = pt + iqz, so that 

Z = * + i <S>= 1+pt + tqZ (13.10) 
1 — pt — iqz 

where p 2 + q 2 = 1. According to (12.36), the first metric function is then 
given by 

x = ^i-tf^l-z 2 

and, from (12.34), u> is given by 

«. = *^£U, W( =(l^!l* s (13.12) 



where 



2az 

$ = - £ — . (13.13) 

(l-pt) 2 + q 2 z 2 V ' 



With this, (13.12) becomes 



_ 4pq(l-t 2 )(l-pt)z _ 2q(l-z 2 )((l-pt) 2 -q 2 z 2 ) 

Uz ' (l-p 2 t 2 -q 2 z 2 ) 2 ' Ut {\-pH 2 -q 2 z 2 ) 2 

(13.14) 

which can immediately be integrated to give 

2 ( p(\ - t 2 )(\ - pt) \ 



where c is an arbitrary constant. 

Chandrasekhar and Xanthopoulos have chosen the constant of inte- 
gration in (13.15) such that uj = when t = 1. The reason for this choice 
is associated with the convenience of aligning the coordinate directions 
with the shear axes on the surface / + g = on which the two waves 
mutually focus each other. Accordingly we set 



c= -1. 



(13.16) 
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A variation of this constant is equivalent to making an Ehlers transfor- 
mation as described in (12.16-18). With (13.16), u may now be written 
in either of the forms 

2q (1 -t)(l -pt - z 2 -pz 2 ) 

ui — — 



(1 + p) (l—p 2 t 2 — q 2 z 2 
2q 2q(l - z 2 )(l - pt) 



2w, (13.17) 



p(l + p) p(l — p 2 t 2 — q 2 z 2 ) 
For purposes of simplification, it is convenient to introduce the terms 
X = (l-pt) 2 + q 2 z 2 

v i 2.2 2 2 2/i .2\ i 2/i 2\ (13.18) 

Y = 1 — p t — q z = p (1 — t ) + q (1 — z ). 
In view of (12.43) and (13.1), in this case we must have 

e~ u = y/l - t 2 \A ~ z 2 = 1 - u 2 - v 2 . (13.19) 

The metric function M can also be obtained from (12.42) and (13.1), 
giving 

e~ M = - X - =. (13.20) 

V 1 — u 2 yl — v 2 

These expressions now complete all the metric functions contained 
in the line element (11.4), which may now be written as 

2Xdu dv Y ( . 2q 2 



ds = — — — — — — I — 

VT^VT^ x\ p(i+p) 



Aq(l-z 2 )( l-pt) ( A _ 2q 
pX 

p 2 XY 



( d * - wh A ») d " <13 - 21) 

(p 2 (l - t 2 )X 2 + 4, 2 (1 - 2 2 )(1 - p<) 2 )d!, 2 . 



After some rearrangement, this may be written in the more convenient 
form 

2 X ( dt 2 dz 2 \ Y 2q x 2 



ds = ^\T^T 2 -T-T 2 ) ~ x [ dx ~ pJiTp) dy 

Ml-^(1-Pt)f dx _2q d \ dy (13 . 22) 



px V p(i + p) 

(1-z 2 ) 



( ((1 - pt) 2 +q 2 ) 2 + p 2 q 2 (l - t 2 )(l - z 2 ))dy 
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in which Y does not appear in a denominator. 

In order to show that this is indeed the Kerr solution, we first relabel 
the parameters by putting 

p = -^M 2 -a 2 /M, q = a/M (13.23) 

which satisfies the condition p 2 + q 2 = 1 for arbitrary constants M and 
a. It is also appropriate to introduce the new coordinates r and 9 instead 
of t and z given by 

t=(r- M)/^/M 2 -a 2 , z = cos9 (13.24) 

where it may be noticed that, using the previous notation, A = tt/2 — 9. 
The remaining coordinates may also be transformed by putting 

T= -^ u ^7wh v - (13 ' 25) 

It is also convenient to introduce the definitions 

p 2 = r 2 + a 2 cos 2 9, A = r 2 -2Mr + a 2 . (13.26) 

With these substitutions the line element (13.22) may then be written in 
the form 

2M 2 ds 2 = 11 2~ J dr 2 3— sin 2 0dTd(/> 

(13 27) 

— I r + a 2^ — Ism 9dq> — p I — dr + dr J 

which after rescaling ds is the standard form of the Kerr solution. 

For this solution to describe the interaction region of colliding waves, 
the coordinates must satisfy the inequality \z\ < t < 1. With (13.24), this 
implies that 

-(M 2 - a 2 ) sin 2 # < A < (13.28) 

which is satisfied only by the region of the Kerr space-time that is inside 
the ergo-sphere. This is not unexpected, as this is the only region of Kerr 
space-time in which the Killing vectors are both space-like. 

The plane on which the two waves collide is given by t = 0, z = 0, or 
by r = M, 9 = n/2 which is between the two horizons. The hyper surface 
on which they mutually focus each other, which is given by / + g = 
or t = 1, is here given by A = which determines the horizons of the 



13.3 The Chandrasekhar-Xanthopoulos solution 



125 



Kerr solution. Thus, the curvature singularity which usually occurs when 
/ + g = 0, in this case, is replaced by a coordinate singularity which is 
normally interpreted as a horizon of the Kerr solution. 

The components of the Weyl tensor inside region IV can be derived 
using (11.10) and have been evaluated by Chandrasekhar and Xanthopou- 
los. The space-time is of algebraic type D and, in particular, it may be 
noted that 



The curvature is clearly non-singular as t — > 1 except in the aligned case 
p = 1, q = (or a = 0). In this case only, the surface t = 1 corresponds 
to a curvature singularity. 

In the interaction region t is a future pointing time-like coordinate. 
It thus follows from (13.24) that r is also a future pointing time-like 
coordinate which increases from the collision at r = M to the outer 
horizon at r = M + \/M 2 — a 2 . In this case a possible further extension 
beyond the horizon would present no problems. 

It may be noted, however, that if the signs of the expressions for t, p 
and q in (13.24) and (13.23) are changed, the Kerr metric (13.27) is still 
obtained but the orientation of r is altered. In this case, the collision at 
r = M is followed by the inner horizon at r = M — \/M 2 — a 2 and any 
extension would lead to a future time-like curvature singularity at r = 0. 

The difference between these two solution has been further analysed 
by Hoenselaers and Ernst (1990). It is the latter case that has been 
described by Chandrasekhar and Xanthopoulos (19866) in which there 
is an extension beyond the horizon at t = 1 which contains a time-like 
singularity analogous to the ring singularity 1 of the Kerr solution at r = 0, 
9 = 7r/2. Since this singularity is time-like, it would be missed by most 
observers. In the other case, which has been given explicitly above, there 
is an analytic extension which is the asymptotically flat exterior Kerr 
solution. It may be pointed out, however, that in both cases the extension 
beyond the horizon is not unique. 

When considering the question of the character of the coordinate 
singularity at t = 1, Chandrasekhar and Xanthopoulos (19866) have 
given a transformation by which it can be removed. They have found 
it convenient to introduce a decreasing time-like coordinate given by 
s = 1 — u 2 — v 2 . The transformation 




(13.29) 



* 2 = 



(1 — pt — iqz) 3 (r — iacos9) 3 



£ = s e 



+ x/q 



C = se 



x/q 



(13.30) 



For a description of the structure of the Kerr solution, see Hawking and 
Ellis (1973), or Chandrasekhar (1983). 
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then removes the coordinate singularity at s = 0, or t = 1. This provides 
an analytic extension beyond the Cauchy horizon which is another part 
of Kerr space-time. This extension, in one case, is followed by a time-like 
curvature singularity which corresponds to the source of the Kerr solution. 

It may be seen that the singularity structure of this solution is very 
similar to that of the degenerate Ferrari-Ibahez solution described in de- 
tail in Section 10.5. In fact this solution reduces to either of the two 
cases of the degenerate Ferrari-Ibanez solution in the limit as q — > 0, 
which is the limit as the polarization of the approaching waves become 
aligned. The non-singular case occurs when p = — 1. It may be noted 
that in this limit a — > 0, the rotation of the Kerr solution vanishes, and 
the solution reduces to the same part of the Schwarzschild space-time as 
the degenerate Ferrari-Ibanez solution. 

It is now appropriate to consider the extension of the solution into 
the prior regions I, II and III which describe the approaching waves that 
give rise to this particular interaction. The metric (13.21) may easily be 
extended into region II simply by replacing both t and z by u. It turns 
out to be convenient to make the coordinate transformation 

x = x+-^y. (13.31) 
1+p 

With this, the line-element in region II takes the form 

IX ( 1 \ 

ds 2 = - du dv - (1 - u 2 ) Xdy 2 + — (dx - 2qu dy) 2 (13.32) 

vi-« 2 V x j 

where X = 1 — 2pu + u 2 . The line-element for region III is identical 
to (13.32) except that v replaces u. Both these line-elements are then 
continuous with the Minkowski space (3.6) in region I. 

The plane wave metric given by (13.32) has the single component of 
the curvature tensor given by 

*4 = -(p-» g )J(tQ- 3 ^ X ~ 2i ^ {1 - ( PU -.\ q f ®{u). (13.33) 

It may immediately be observed that the approaching waves have variable 
polarization except in the aligned limit q — > 0. They include an impulsive 
component and a step component. 

As in the degenerate Ferrari-Ibanez solution, the singularity struc- 
ture is thus complicated by the fact that the solution contains impulsive 
waves on the wave fronts u = and v = which form the boundaries of 
the interaction region. The presence of these wave components introduces 
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additional distribution valued singularities at the points u = 0, v = 1 and 
v = 0, u = 1. As described in Chapter 8, these point singularities ensure 
the existence of fold singularities on the surfaces u = 1 and v = 1 in 
regions II and III respectively. 

The profiles for the approaching waves in both cases have been de- 
scribed in more detail by Hoenselaers and Ernst (1990). They have shown 
that, in the case where the Cauchy horizon corresponds to the inner Kerr 
horizon, the amplitude diverges towards the fold singularity which is thus 
a non-scalar curvature singularity. On the other hand, when the Cauchy 
horizon corresponds to the outer Kerr horizon, the amplitude decays to- 
wards a quasiregular fold singularity. 

13.4 Other solutions 

It may be noticed that the solutions presented in the previous two sec- 
tions have been obtained using the generating techniques described in 
Sections 12.3 and 12.4 respectively. Other techniques may also be used, 
both in isolation and in various combinations. Since they involve the same 
Ernst equation, it is clear that all the solution generating techniques that 
have been developed for stationary axisymmetric space-times in recent 
years can also be applied to colliding plane wave situations. 

A number of papers have recently been published which present new 
colliding wave solutions using some of these techniques, and it is likely 
that many more such papers will appear. Most of these concentrate on a 
description of a particular technique and how it may be adapted to the 
colliding wave situation. At present, however, there is less emphasis on the 
physical significance of the particular solution generated. Nevertheless, 
solutions with a number of parameters have been obtained, and general 
classes of solutions have been defined. 

Ideally, one would like to be able to specify the particular solution 
corresponding to an arbitrary set of initial conditions. Such techniques 
are available in the axisymmetric case, and may soon be extended to 
the colliding wave situation. The initial attempts at the solution of this 
problem are mentioned in the next chapter. 

So far most of the explicit solutions that have been presented have 
involved impulsive wave components. In the notation of (7.11), they have 
used ni = 77-2 = 2. Much less attention has been given to smooth-fronted 
waves of finite duration. 

Some further general aspects of cases involving impulsive components 
have been clarified by Ferrari (1988). In particular, it has been shown that 
the proper time between the collision and the subsequent singularity is 
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inversely proportional to the square root of the amplitudes of the impul- 
sive components of the approaching waves. This time is also affected by 
their relative polarization and is a minimum when the approaching waves 
are colinear. This result also applies to solutions which contain a Cauchy 
horizon rather than a curvature singularity. In addition, Ferrari has also 
shown that the shock waves which accompany the impulsive waves only 
affect the rate at which the Weyl scalars diverge on the singularity. 

Having made these general remarks, we may now briefly review the 
new exact solutions that have been presented. The generation techniques 
themselves will not be described in detail here. 

First, there is the important class of solutions that has been ob- 
tained by Ferrari, Ibahez and Bruni (1987a, b) using the inverse scatter- 
ing method developed by Belinskii and Zakharov (1978, 1979) and Carr 
and Verdaguer (1983). Using this method, they have obtained a two- 
parameter class of colliding wave solutions with non-aligned polarization. 
This in fact is a generalization of the soliton solution of Ferrari and Ibanez 
(19876) described in Sections 10.4 and 10.5. 

The Ferrari-Ibanez-Bruni solution has been obtained using as seed 
the Kasner metric. They utilize the Belinskii-Zakharov soliton technique 
with two real poles. The resulting metric for region IV is singular on 
the space-like surface / + g = 0, except in a particular limiting case. By 
extending the solution back into regions II and III, it can be seen that the 
approaching waves must have both an impulsive component and a step 
component with variable polarization. 

The limiting case of the Ferrari-Ibanez-Bruni solution is particularly 
interesting. It is of type D, and can be shown to be part of the Taub- 
NUT solution 2 in the Taub region where there are two space-like Killing 
vectors. The line element can be written in the form 

ds 2 = C(l + 2psinV> + sin V> 2 )(d-i/> 2 - dA 2 ) 

- *T Si °^. 2 , ) (da: - 2qsm\dy) 2 (13.34) 
V 1 + 2p sin ip + sin ip J 

— cos 2 A(l + 2psin^ + sin 2 ip)dy 2 

where p and q are constants satisfying p 2 + q 2 = 1. It can immediately be 
seen that, when q = and p = ±1, this further reduces to the two colinear 
degenerate Ferrari-Ibahez solutions that have already been considered in 
Section 10.5. The structure of this solution has been further examined by 
Ferrari and Ibanez (1988). 

2 Newman, Tamburino and Unti, (1963); see also Hawking and Ellis (1973) 
section 5.8. 
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Ernst, Garcia-Diaz and Hauser (1987a) have shown that an Ehlers 
transformation (12.13) applied to the Ernst function of the Nutku-Halil 
solution (13.3), together with a simple coordinate transformation, yields 
a non-colinear generalization of the Ferrari-Ibanez soliton solution that 
was described in Section 10.4 but only with even integer values of the 
parameter a. This generalization is in fact contained in the solution of 
Ferrari, Ibahez and Bruni (1987a, b). 

A further generalization extending the Ferrari-Ibanez-Bruni solu- 
tion was reported by Ernst, Garcia-Diaz and Hauser (1987a, 6), with a 
full derivation given in Ernst, Garcia-Diaz and Hauser (1988). This so- 
lution contains three arbitrary parameters. In stating this solution, it is 
convenient first to define the function 



1 (a, a, p) = 



(p+p ">(Wt) +(p - p,) (r^) ' 



+ 



i\J\ — z 2 



\ «/2 /i , \ a/2' 



(13.35) 



where a, a and (3 are the arbitrary parameters, and where 

p = cosa, q = sin a, p' = cos/3, q' = sin(3. (13.36) 

With this definition, the Ernst function for this new solution may be 
written as 

Z(a,a,P) = (1 - t 2 r /2 (l ~ z 2 r /2 l { ( a + y« a \ . (13-37) 

1 {a — 1, p, a) 

When a = 0, this reduces to the Nutku-Halil solution. With \a\ = 1, 
one obtains the Kerr-NUT solution with a being the NUT parameter. 
The subcase with \a\ = 1 and f3 = becomes the Schwarzschild-NUT 
space-time. 

Transforming (13.37) using the Neugebauer-Kramer involution pro- 
duces another three-parameter family of solutions with the alternative 
Ernst function given by 

Z(a,a,0) = (1 - t 2 )^/ 2 (l - z 2 ) (1 -° )/2 r( °~ 2,a '^ . (13.38) 

T(a, a, p) 

With a = 0, this reduces to the Chandrasekhar-Xanthopoulos solution. 

Ernst, Garcia-Diaz and Hauser (19876) have also found a simpler 
realization G (Xi) of the Geroch group of transformations relating collid- 
ing wave solutions. With this , they have obtained a generalization of the 
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above solution in the case when a = 2. A further generalization of the 
a = 3 solution is also obtained using a Neugebauer-Kramer involution. 

It was pointed out by Ernst (1986) that the Nutku-Halil solution 
can be derived from the isotropic (n = 0) Kasner metric by apply- 
ing a double-Harrison (Backhand) transformation. The Kerr space-time 
of the Chandrasekhar-Xanthopoulos solution can also be derived from 
Minkowski space using a double-Harrison transformation. It therefore 
seemed reasonable that other colliding wave solutions may be generated 
by applying double-Harrison transformations to other Kasner metrics. 
This procedure has indeed been successful. It was used by Ernst, Garcia- 
Diaz and Hauser (1988) to determine the solution described above in 
(13.35) to (13.38). 

The effect of a double-Harrison transformation can be deduced by 
any number of competing techniques. Naturally, however, Ernst, Garcia- 
Dfaz and Hauser preferred to use the so-called Hauser-Ernst homogeneous 
Hilbert problem approach (see also Ernst 1988). 

Finally, it may be mentioned that Hassan, Feinstein and Manko 
(1990) have recently obtained another particular solution describing the 
non-colinear case using the generating technique of Gutsunaev and Manko 
(1988) as described here in (12.56) and (15.57). Further solutions could 
easily be obtained using this method. 



